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Interior Radiances in Optically Deep Absorbing Media 
II, Rayleigh Scattering 

Gilbert N, Plass, George W« Kattawar, and Judith Binstock 
Department of Physics, Texas A&M University 
College Station, Texas 77843 

ABSTRACT 

The interior radiances are calculated within an optically deep ab- 
sorbing medium scattering according to the Rayleigh phase function. The 
accuracy of the matrix operator method is improved by many orders of mag- 
nitude through the use of accurate starting values obtained by the Runge- 
Kutta method rather than from the single scattering approximation. The 
radiance and flux are given for a range of solar zenith angles and for 
single scattering albedos of 1, 0.99, 0.9, 0.5, and 0.1. The development 
of the asymptotic angular distribution of the radiance is illustrated. 

It is shown that this asymptotic distribution is probably physically unob- 

_8 

servable when w o <0.8, since the flux is less than 10 of its original 
value at the beginning of the asymptotic region. The ratio of the upward 
to downward flux is calculated and is shown to be remarkably constant with- 
in the medium except very close to the boundaries. The heating rate with- 
in the medium is found to be very nearly proportional to the downward flux, 
except near the boundaries. When the single scattering albedo is small, 
a number of examples illustrate the significant contribution of the direct 
solar flux to the total flux even at great optical depths within the medium. 
The total downward flux decreases exponentially with optical depth away 
from boundaries when the single scattering albedo is greater than or equal 
to 0.9; when it is less than or equal to 0.5 only an approximate exponential 
fit can be obtained within the region accessible to experimental observa- 
tion. 
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INTRODUCTION 

Results are given in this paper for the interior radiance within 
optically deep absorbing clouds and its dependence on the single scatter- 
ing albedo and the optical depth within the cloud. The matrix operator 
theory recently reviewed by Plass, Kattawar, and Catchings^ is used to 
obtain an entirely rigorous numerical solution of the radiative transfer 
equations. The Rayleigh phase function is used for illustration in these 
first calculations, although the method is applicable to general phase 

functions. Further details of the matrix operator method are given in the 

( 2 ) 

first part of this paper by Kattawar and Plass (hereafter referred to 
as I). A starting proceedure which uses the Runge-Kutta method instead 
of the usual single scattering approximation improves the accuracy by 
many orders of magnitude in typical cases as described in I. A detailed 
discussion of errors in a one-dimensional matrix operator calculation is 
presented in I; in most cases these results apply equally well to the 
multidirectional case. 

Interior radiances can easily be calculated by Monte Carlo methods. 

(3 

Results obtained in this manner have been reported by Plass and Kattawar 
and Kattawar and Plass Interior radiances were obtained at points 

within the earth's ocean-atmosphere system. These appear to be the only 
reasonably accurate values for the interior radiance which have been re- 
ported in the literature. A relatively simple equation for the interior 
radiances originally given by Bellman ^ (see equations (3-4) on p. 348) 
and based on the matrix operator theory does not appear to have been used 
previously in practical calculations. 

The results given here show the dependence of the interior radiance 
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and flux on the optical depth within the medium as well as on its absorb- 
ing and scattering properties. The development of the asymtotic angular 
distribution of the radiance is illustrated. This asymptotic distribution 
is unobservable when the single scattering albedo is small, since the 
region only begins at great optical depths where the flux is extremely 
small. The range of optical depth for which the ratio of the upward to 
downward flux is constant is investiaged as is the region in which the de- 
crease of the downward flux with optical depth can be represented by an 
exponential. The heating rate within the medium is calculated and found 
to be nearly proportional to the downward flux except near boundaries. 

DOWNWARD RADIANCE 

The interior radiances were calculated from Equations (5) and (6) of 

Plass, Kattawar, and Catchings^ by the methods of matrix operator theory 

as reviewed by them. Greatly improved accuracy was obtained by the use 

of a Runge-Kutta method to obtain a starting value instead of using the 

single scattering approximation. The differential equations satisfied by 

the reflection and transmission operators (see Kattawar were integrated 

-3 

from the origin to an optical depth of the order of 10 . Since the error 

of the Runge-Kutta method is proportional to h^ where h is the interval 
size, the solution at this optical depth has an error of the order of 10 
This is approximately equivalent to machine accuracy since all calculations 
reported here were done in double precision on the IBM 360/65 (the equivalent 
of 16 significant figures). Each of the five cases reported here required 
approximately 4 minutes of computer time. 

The calculations were continued out to very large (x = 16,794) opti- 
cal depths when the single scattering albedo = 1; when a> o < 1 the 



4. 


calculations were carried out to optical depths of 45 to 109 depending on 
the particular value of In each case results were calculated for 

several different values for the albedo of the lower surface. However, 
this is not of particular interest for this problem as all of the results 
presented here are for points sufficiently far away from the lower bound- 
ary so that the radiance values are independent of the surface albedo. 

The downward normalized diffuse radiance is shown in Fig. 1 when 

to =0.99 and to = 0.9 as a function of the cosine y of the zenith angle 
o o 

of observation. These curves are in the principal plane, so that the 
azimuthal angle <f> = 0° or 180°. The cosine y Q of the solar zenith angle 
is 0.85332 which corresponds to a zenith angle 0 = 31.42°. The downward 
radiance is shown for a number of different values of the optical depth 
x within the medium. In each case the radiance is multiplied by the fac- 
tor (ir/diffuse flux at depth t), so that variation of the radiance with y 
at different depths can conveniently be compared. The incoming flux is 
normalized to unity across a plane at right angles to the incoming beam. 

When id =0.99, the downward radiance is a maximum at the horizon 
o 

and has a minimum value just beyond the zenith on the antisolar side when 
t < 1. When x = 8.9795, the radiance has reached its asymptotic form 

(Q\ 

predicted by Preisendorfer v ' . The curves for higher values of t are not 

plotted, since they are identical with this last curve to the accuracy 

that can be shown in the figure. The asympotic radiance is, of course, 

symmetric around the zenith with a maximum value at the zenith. 

The downward normalized interior radiance for w = 0.9, 0.5, and 0.1 

o 

is shown in Figs. 1, 2, and 3 respectively. The angular dependence of 
the radiance when x << 1 changes very little as u) decreases. On the 

other hand the angular dependence for large values of x depends critically 
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on a) . The ratio of the radiance at the zenith to that at the horizon 
o 

increases greatly as us q becomes smaller. The asymptotic form for the down- 
ward radiance is valid approximately when x > 9 provided that w o >_ 0.9. 

On the other hand when ur ^0.5 the asymptotic form has not developed even 
at t = 40. Although the existence of an asymptotic form has been predicted 

(Q\ 

by Preisendorfer v , the optical depth required for this angular distri- 
bution has not been known previously. 

The variation of the flux and radiance as a function of optical depth 
is given for each case in Tables 1-6. Most of the quantities in these 
tables are discussed in a later section. The optical depth is given in the 
first column. In each case the lower boundary of the medium was taken at 
the largest optical depth shown in each table. A lower boundary surface 
of zero albedo was assumed. The seventh column • ' •• of Tables 1 and 6 

and the eighth column of Tables 2-5 give the ratio of the downward dif- 
fuse radiance at the zenith (y = 1) to that near the horizon (p = 0.037850 
or 0 = 87.83°) when <J> = 0° and p q = 0.85332 (31-42°). When = 1 (Table 
1), this ratio is small near the upper boundary, becomes very slightly 
larger than unity in the interior of the medium, and then further increases 
to 2.71 at the lower boundary. When o) q = 0.99 (Table 2), the ratio has 
the limiting value of 1.20451 which is accurate to six significant figures 
when 16 < x < 69. When cl) q = 0.9 (Table 3), the asymptotic value is 

2.13053 which is calculated to six significant figures when 32 < x < 65; 

the ratio is within 10% of the asymptotic value when 8 < x < 73. 

However, the situation is quite different when u) q = 0.5 and 0.1 

(Tables 4 and 5). The ratio continues to increase as x increases in each 
of these cases and no asymptotic angular form develops when x < 45. Of 
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course, an asymptotic form would develop at sufficiently large values. 

However, it should be pointed out that the total downward flux is only 
-19 -20 

7.5 x 10 and 3.8 x 10 for m = 0.5 and 0,1 respectively at x = 45. 

The flux is normalized in these calculations to unit incoming solar flux 

per unit area perpendicular to the incoming beam. Thus when ^ <0.5, 

o — 

the asymptotic angular form has not developed at optical depths such that 

-20 

the downward flux has been reduced to 10 of its original value. Ex- 
perimental observation of the asymptotic form is out of the question under 
these conditions. 

All of the values presented in Tables 1-5 assume that the sun is 

at the zenith. Selected values for other solar zenith angles are given in 

Table 6. The ratio of the downward flux already discussed is the same to 

six significant figures in the asymptotic region for y Q = 0.18816 

(0 “79.15°) as for y = 1.0, when u> “1, 0.99, 0.9. On the other hand 

o o 

when 0 )^ = 0.5, the ratio depends significantly on the solar zenith angle 

out to optical depths of the order of 50, as is shown in Fig. A. 

Even when x = 50 the ratio has a value that is an order of magnitude 

smaller when y = 0.18816 than when y = 1 for the case w =0.1. From 
o o o 

the slow convergence of these curves it would appear that an optical depth 

of the order of one thousand would be required before these curves would 

approach a limit and the asymptotic region would be reached. The flux 

-435 

would be of the order of 10 at such large optical depths and completely 
unobservable. 

( 8 ) 

The prediction of Preisendorfer v ' that the asymptotic angular dis- 
tribution becomes highly peaked in the forward direction as the single 
scattering albedo decreases is confirmed. However, our calculations show 
that the depth at which the asymptotic distribution develops also increases 


7. 


as the single scattering albedo decreases. Even for the case cd o = 0.9, 

“8 

the downward flux is 4.5 x 10 of its original value at the beginning of 
the asymptotic region. Radiance values of this order cannot be detected 
experimentally against the background of other radiation sources. Thus, 
there is some limiting value of u) q such that it becomes exceedingly diffi- 
cult to observe the asymptotic radiance distribution for smaller values of 
w . From the calculations given here it might be presumed that the limit- 
ing value is approximately = 0.8. 

All of the radiance values presented so far have been in the principal 
plane which contains the incident solar direction. Some examples of the 
azimuthal variation of the downward radiance are given in Figs. 5 and 6. 

As before the radiance has been multiplied in each case by the factor 
(tt / diffuse flux at depth x), so that the variation of the radiance with u 
at different depths can conveniently be compared. The values p q = 0.85332 
(31.42°) and p => 0.71392 (44.44°) were chosen, since a fairly large azi- 
muthal variation occurs for these values. When w =0.5 the downward 

o 

normalized interior radiance decreases as 4> increases when t is small. It 
reaches a maximum value when t is about 3 , and exhibits less and less de- 
pendence on <|> as t increases further. When t = 32, no variation with azi- 
muth can be seen on the scale of the figure. The variation is qualitatively 
similar when u) q = 0.1, except that the variations are more pronounced. The 
curve is still decreasing when x = 41, since the asymptotic region has not 
yet been reached. 

UPWARD RADIANCE 


The upward normalized interior radiance is shown in Figs. 7, 8, and 9 
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when a> o = 0.99, 0.5, and 0.1 respectively. The angular variation of the 
upward radiance is very much less than that of the downward and has been 
plotted on a linear scale. The radiance is these figures is shown as a 
function of the cosine y of the nadir angle and is for the incident plane 
(<(> = 0° or 180°). The solar horizon is on the left of all the figures. 

The cosine of the solar zenith angle was chosen as = 0.85332 (31.42°). 

In all cases the upward radiance has been normalized by multiplacation by 
the factor (ir/diffuse flux at depth t), as was previously explained. 

When = 0.99 and the total optical depth of the medium is 109, the 
upward radiance near the upper surface of the medium has a minimum at the 
solar horizon, increases through the nadir, reaches a maximum value near a 
nadir angle of 50°, and then decreases toward the antisolar horizon. The 
upward normalized radiance increases at the horizon and decreases toward the 
nadir at greater optical depths within the medium. A limiting angular de- 
pendence of the normalized radiance is reached at about t = 9 with a mini- 
mum at the nadir and a maxima at both the solar and antisolar horizons. 

When oj q = 0.5, Fig. 8 shows that the normalized upward radiance has a 
somewhat different form at small optical depths. It decreases from the 
solar horizon to a minimum value at a nadir angle of about 40° and then 
decreases uniformly through the nadir to the antisolar horizon. As the 
optical depth increases, the normalized upward radiance undergoes only a 
small change in its angular variation into its limiting angular dependence 
which must be symmetric around the nadir. This limiting form occurs at 
about t = 25. 

The normalized upward radiance when = 0.1 is shown in Fig. 9. 

The angular variation at small optical depths has the same qualitative 
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variation as was shown for the previous case. The asymptotic angular var- 
iation has not been reached even at t = 42, although the curve is beginn- 
ing to show only a small variation with t. In this case the asymptotic 
curve appears to have a minimum near a nadir angle of 55° and a relative 
maximum at the nadir and the horizon. 

The variation of the normalized upward radiance with azimuthal angle 
is relatively small as is shown in Fig. 10. The following parameters were 
selected for this illustration: y Q = 0.85332 (31.42°), y = 0.71392 

(44.44°), and u = 0.5. At small optical depths within the medium the 
radiance increases slowly as <p increases from 0° and 180°. There is no 
appreciable azimuthal variation when t > 17 . 

The sixth column of Tables 1 and 6 and the seventh column of Tables 

2-5 give the ratio of the upward radiance at the nadir (y = 1) to that 

near the horizon (y = 0.03785 or 0 8 87.83°). The limiting value for this 

ratio holds over approximately the same range of optical depths as for the 

corresponding ratio for the downward radiance already discussed. The 

limiting value for the ratio for the upward radiance in the interior of 

the medium is 0.99994, 0.86137, and 0.69386 when u) = 1.0, 0.99, and 0.9 

o 

respectively. The ratio for the upward radiance does not vary over as wide 
a range as that for the downward radiance. For example, when w o = 0.1, 
the ratio for the upward radiances is 1.0300 at the upper boundary and has 
changed only to 0.9422 at an optical depth of 40.98, while the correspond- 
ing ratio for the downward radiance varies from 0.075 to 75.6. 

FLUX 

Flux values at various optical depths within the medium as given in 
Table lwhen u ° 1. The second column gives the upward flux; the third 
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column shows the downward diffuse flux; the fourth column gives the total 
downward flux including the incident solar beam which is assumed to be at 
the zenith. The difference between the total downward flux and the upward 
flux is shown in the fifth column. The upward and downward flux have the 
same numerical value to four or five significant figures. Nevertheless 
their difference is constant to four significant figures in our calcula- 
tions. The difference is necessarily a constant for a conservative problem 

(w =1). The incident solar beam makes no contribution to the total down- 

o 

ward flux to five significant figures at optical depths greater than 13. 

The total downward flux becomes significantly greater than unity over a 

large range of optical depths starting from the upper surface, a phenomenon 

which does not violate the conservation of energy. The explanation is 

(9) 

similar to that given by Plass and Kattawar for the fact that the down- 
ward flux from the lower surface of a medium may be greater than unity under 
certain conditions. 

The first four columns in Tables 2-5 tabulate the same quantitites 
as are in Table 1. The fifth column of Tables 2-5 shows the ratio of the 
upward flux to the total downward flux. The quantity 3F/9 t, which is pro- 
portional to the heating rate of the layer, is given in the sixth column. 
Within the tabulated accuracy of four significant figures, the incident 
beam does not contribute to the total flux at optical depths greater than 

17 when <o = 0.99 or 0.9. However, the situation is quite different when 
o 

-9 

co =0.5. The total downward flux is 7.5 x 10 at t = 45 and yet the 
o 

direct beam is still contributing 4% of the total flux value. Similarly 

-20 

when (o =0.1, the total downward flux is 3.4 x 10 at t = 45, so that 
o 

the direct beam is contributing 77% of the total downward flux. Since the 
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contribution of the direct beam to the total downward flux is 2.9 x 10 
at this optical depth, this is a lower limit for the total downward flux no 
matter how small may be. When 0.5 and y = 1, a significant number 

of the photons occurring at great depths have been a part of the unscattered 
direct beam until these depths are reached. 

The ratio of the upward interior flux to the total downward interior 
flux is often measured experimentally. When = 0.99, this ratio is re- 
markably constant except near the lower boundary. The tabulated values 
range only from 0.753 to 0.793. The ratio has the constant value 0.793 when 
9 £ t _< 85. Similarly the ratio has the constant value 0.470 when 
17 £ t £ 65 and = 0.9. This ratio only varies over a range of 2%, 5%, 

5%, and 1% from an optical depth of unity to near the lower boundary of the 
medium for = 0.99, 0.9, 0.5, 0.1 respectively. The physical explanation 
for the constancy of this ratio is that the upward flux is derived by mul- 
tiple scattering from the downward flux, so that the ratio varies with op- 
tical depth within the medium only as other factors change the relative con- 
tributions of the scattering above and below the point that is being con- 
sidered. Thus the ratio would be expected to be constant in the region where 
the asymptotic radiance distribution holds as well as approximately constant 
at all optical depths except those very close to boundaries. 

The quantity 9F/9x, which is proportional to the heating rate, is also 
given in Tables 2-5. In all cases it is proportional to the total down- 
ward flux when x > 1 and the interior point is not too near the lower bound- 
ary. In all cases the heating rate increases as the optical depth increases 
starting from the upper boundary and reaches a maximum when x is approxi- 
mately 1, 0.25, and 0.01 when = 0.99, 0.9, and 0.5 respectively. 


12 . 


Since all of the data in Tables 1-5 are for the case w =1, some 

o 

selected data for other values of y are shown in Table 6. The ratio of 

o 

the upward to the total downward flux is given in column five. This ratio 
in the asymptotic region is 0.793 and 0.470 for = 0.99 and 0.9 respec- 
tively when p Q = 0.18816 (0 q = 79.15°). This agrees exactly with the 
asymptotic value given in Tables 2 and 3 for y Q = 1. When = 0.5 the 
asymptotic value for this ratio differs only in the fourth significant fig- 
ure for different values of y (Tables 4 and 6). When to =0.1, this ratio 

o o 

had not approached its limiting value even at t =41. The value of the 

ratio at x = 41 is 0.01805, 0.01826, and 0.01820 for y =1, 0.8332, and 

o 

0.18816 when u> = 0.1. In the true asymptotic region the value of this 

ratio should be independent of y Q . 

The ratio of the direct downward flux to the total downward flux is 

given in the fourth column of Table 6. It is important to understand the 

range of optical depths over which the direct flux is an important part of 

the total flux. For example, when a) Q = 0.1 and y Q = 1, Table 5 shows that 

the direct flux is more than 10% of the total for all x < 45. From Table 

6 and additional data from the original calculation we find that the direct 

flux is more than 10% of the total when x < 33 for y = 0.85332 (31.42°) 

o 

and when x < 1.5 for y = 0.18816 (79.15°). Thus in this case when the 

o 

single scattering albedo is small the importance of the direct beam depends 

critically on the solar zenith angle. The exponential degradation of the 

solar beam measured along the vertical is, of course, much more rapid when 

the sun is near the horizon increasing as sec 6. It is also interesting 

to study how rapidly thi9 ratio decreases with optical depth for other 

values of to from Table 6. 
o 
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The total downward flux decreases exponentially when the optical depth 

is sufficiently great and the point of observation is not too close to the 

lower boundary. The decrease with optical depth of the total downward flux 

can be represented quite accurately with exp (-b-r), where b = 0.17244 and 

0.52318 for u) q = 0,99 and 0.9 respectively, provided that t > 3 and is not 

near the lower boundary. However, when the absorption is relatively large 

(w £ 0.5), only an approximate fit can be obtained with b = 0.944 when 

to = 0.5 and b = 0.997 when u> = 0.1. The reason for the approximate fit 
o o 

is two-fold: (1) the direct beam is a significant fraction of the total 
downward flux even up to t = 40 for normal incidence with large absorption. 
Since the direct beam decays at a different rate than the diffuse flux, the 
latter rate cannot be established until the direct beam is no longer con- 
tributing a majority of the photons in the vicinity. (2) When the absorp- 
tion is large, a majority of the photons at great optical depths arise from 
single scattering events at smaller optical depths which scatter the photon 
out of the direct beam into a nearly vertical direction. The number of 
photons which have undergone single scattering and reached large optical 
depths decreases with optical depth as exp(-j). Deviations from this rate 
of decrease for the total flux are due to photons that have undergone higher 
orders of scattering. An equilibrium between the number of photons that 
have undergone various orders of scattering is not obtained in the case of 
large absorption until very large optical depths are reached. 

The development of the exponential decrease in both the flux and rad- 
iance is shown in Fig. 11 for the case when w =0.1. The radiance or flux 

o 

is shown as a function of the optical depth within the medium. In each case 
the dashed line represents a pure exponential dependence obtained by fitting 
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the curve at large optical depths. On the rather small scale of this fig- 
ure deviations are observed up to optical depths of the order of 20. The 
actual calculations show that there are still deviations from the pure ex- 
ponential as large as 1% at t = 40. 

The upper curve in Fig. 11 shows the flux for v = 1, while the next 
curve gives the flux for y Q = 0.18816 (79.15°). In the first case the 
actual flux for t < 20 is below, while in the second case it is above, the 
value given by the pure exponential indicated by dashed lines in the figure. 
The downward radiance is shown for five different values of y (1, 0.92226; 
0.81141; 0.47254; 0.037850) when y Q = 0.18816. The downward radiance near 
the zenith is much larger than that nearer the horizon for large t due to 
the large maximum that develops at the zenith (see Fig. 3). The downward 
radiance near the zenith approaches the limiting exponential curve from 
below as t increases, while it approaches it from above for larger zenith 
angles. For this particular case, the downward radiance at the zenith 
(y = 1) decreases exponentially with the optical depth within the accuracy 
of Fig. 11 when t > 2, while the downward radiance for y = 0.47154 decreases 
exponentially only when t > 20. Thus the pure exponential decrease develops 
at much smaller values of the optical depth for the radiance near the zenith 
than near the horizon. 

CONCLUSIONS 

The interior radiance and flux have been studied for a very thick 
Rayleigh scattering absorbing layer. Their dependence on the solar zenith 
angle, the single scattering albedo, and the optical depth within the medium 
has been shown. The development of an asymptotic angular form for the down- 
ward radiance as predicted by Preisendorfer is confirmed. A range of optical 
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depths is given for each case for which this asymptotic form is valid. 

When < 0.8, the asymptotic form only develops at great optical depths 

—8 

such that the downward flux is 10 or less of its value at the upper 

boundary. Thus in this case the flux is so small that the asymptotic form 

could not be observed experimentally. As an example, when = 0.1, the 

asymptotic form has not developed even at t =41, at which depth the down- 

*“18 

ward flux is 2 x 10 times its value at the upper boundary. 

The ratio of the upward to the downward flux is remarkably constant 
within the medium. This ratio never varies by more than 5% from an optical 
depth of unity to near the lower boundary of the medium for any value of the 
single scattering albedo. Over the same range of optical depths within the 
medium, the heating rate is found to be proportional to the total downward 
flux. 

When >_ 0.9 the total downward flux decreases exponentially with 
optical depth provided x > 3 and is not near the lower boundary. However 
when « _< 0.5, only an approximate exponential fit can be obtained at least 

down to optical depths of 40 where the flux is so small as to be experi- 
mentally unobservable. When _< 0.5, the direct flux makes an important 
contribution to the total flux down to optical depths of the order of 45. 

For example, when m =0.1 and P =1, the direct flux is more than 10% 

o o 

of the total flux for all x < 45. 
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TABLE 1 


■ 0) = 1 , u =1 

o o 

Optical F up F d °™ F d °“", pdown _ F ap up / up I d °™( 1 )/I da > .03785: 

Depth diffuse total total 


er 

0.9999 

0 


1.0000 

1.004 

-4 

1.57779 


0.125 

1.0443 

0.1619 


1.0444 

1.004 

-4 

1.38629 

0.14608 

0.25 

1.0781 

0.2994 


1.0782 

1.004 

-4 

1.29104 

0.24696 

0.5 

1.1282 

0.5217 


1.1283 

1.004 

-4 

1.18434 

0.40630 

1 

1.1888 

0.8211 


1.1889 

1.004 

-4 

1.08845 

0.61795 

5 

1.2635 

1.2569 


1.2636 

1.004 

-4 

1.00096 

0.98839 

13 

1.2639 

1.2640 


1.2640 

1.004 

-4 

0.99994 

1.00005 

77 

1.2591 

1.2592 


1.2592 

1.004 

-4 

0.99994 

1.00006 

205 

1.2495 

1.2496 


1.2496 

1.004 

-4 

0.99994 

1.00006 

8,397 

0.6324 

0.6325 


0.6325 

1.004 

-4 

0.99988 

1.00011 

16,589 

1.544 -2 

1.544 

-2 

1.544 -2 

1.004 

-4 

0.99532 

1.00468 

16,717 

5.802 -3 

5.903 

-3 

5.903 -3 

1.004 

-4 

0.98761 

1.01238 

16,781 

9.824 -4 

1.083 

-3 

1.083 -3 

1.004 

-4 

0.92963 

1.06998 

16,793 

7.820 -5 

1.786 

-4 

1.786 -4 

1.004 

-4 

0.48306 

1.56051 

16,793.875 

1.038 -5 

1.107 

-4 

1.107 -4 

1.004 

-4 

0.12149 

2.29156 

16,794 

0 

1.004 

-4 

1.004 -4 

1.004 

-4 
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0 
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0 
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0.7529 
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-3 
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2.436 

-2 

1.39077 
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-3 

1.0006 
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-2 
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0.05341 

0.01074 

0.7550 
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-2 

1.0021 

0.753 

2.480 

-2 

1.37015 

0.05848 

0.04199 

0.7611 
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-2 
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0.756 

2.592 

-2 
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0.1045 

0.7705 


0.1143 


1.0150 

0.759 

2.746 

-2 
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0.4795 

0.7917 


0.4059 


1.0250 

0.772 

3.161 

-2 

1.05250 

0.41733 

0.9795 

0.7756 


0.6172 


0.9927 

0.781 

3.276 

-2 

0.95979 

0.65973 

2.98 

0.5950 


0.7009 


0.7517 

0.792 

2.652 

-2 

0.87438 

1.07010 

8.98 

0.2145 


0.2703 


0.2704 

0.793 

9.639 

-3 

0.86144 

1.20293 

16.98 

5.399 

-2 

6.807 

-2 

6.807 -2 

0.793 

2.426 

-3 

0.86137 
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24.98 
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-2 

1.713 

-2 

1.713 -2 

0.793 
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-4 
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32.98 
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-3 
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-3 

4.312 -3 
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-4 

0.86137 
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8.611 

-4 
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-3 
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0.793 

3.870 

-5 
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1.087 

-4 

1.371 

-4 
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0.793 

4.887 

-6 
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-6 
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-6 

8.685 -6 

0.793 

3.096 

-7 

0.86137 

1.20451 

84.98 

4.364 

-7 

5.502 

-7 

5.502 -7 

0.793 

1.961 

-8 

0.86130 

1.20457 

92.98 

1.094 

-7 

1.382 

-7 

1.382 -7 

0.792 

4.922 

-9 

0.86028 

1.20557 

.00.98 

2.593 

-8 

3.350 

-8 

3.350 -8 

0.774 

1.181 

-9 

0.84328 

1.22211 

08.98 

0 


3.235 

-9 

3.235 -9 

0.000 

5.560 

-11 
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Depth 

F UP 
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0 

0.4183 

0 


1.0000 

0.418 
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0.000977 

0.4183 

7.984 

-4 

0.9998 

0.418 

0.1822 

1.15990 

0.05681 

0.00293 

0.4183 

2.389 

-3 

0.9995 

0.418 

0.1826 

1.15666 

0.05815 

0.01074 

0.4182 

8.669 

-3 

0.9980 

0.419 

0.1844 

1.14448 

0.06368 

0.04199 

0.4176 

3.262 

-2 

0.9915 

0.421 

0.1898 

1.10484 

0.08844 

0.1045 

0.4148 

7.587 

-2 

0.9766 

0.425 

0.1957 

1.04972 

0.14684 

0.2295 

0.4058 

0.1479 


0.9428 

0.430 

0.2002 

0.97896 

0.27068 

0.4795 

/ 

0.3806 

0.2486 


0.8677 

0.439 

0.1975 

0.89709 

0.48340 

0.9795 

0.3204 

0.3384 


0.7139 

0.449 

0.1749 

0.81625 

0.80526 

2.98 

0.1304 

0.2303 


0.2811 

0.464 

7.553 -2 

0.72391 

1.52815 

8.98 

6.007 -3 

1.267 

-2 

1.280 -2 

0.469 

3.546 -3 

0.69507 

2.07690 

16.98 

9.164 -5 

1.951 

-4 

1.951 -3 

0.470 

5.414 -5 

0.69388 

2.12884 

24.98 

1.394 -6 

2.968 

-6 

2.968 -6 

0.470 

8.237 -7 

0.69386 

2.13048 

32.98 

2.121 -8 

4.516 

-8 

4.516 -8 

0.470 

1.253 -8 

0.69386 

2.13053 

40.98 

3.227 -10 

6.870 

-10 

6.870 -10 

0.470 

1.906 -10 

0.69386 

2.13053 

48.98 

4.909 -12 

1.045 

-11 

1.045 -11 

0.470 

2.900 -12 

0.69386 

2.13053 

56.98 

7.468 -14 

1.590 

-13 

1.590 -13 

0.470 

4.412 -14 

0.69386 

2.13053 

64.98 

1.136 -15 

2.419 

-15 

2.419 -15 

0.470 

6.712 -16 

0.69385 

2.13053 

72.98 

1.704 -17 

3.668 

-17 

3.668 -17 

0.464 

1.014 -17 

0.68433 

2.13978 

76.98 

0 

3.508 

-18 

3.508 -18 

0.000 

5.787 -19 


3.76160 
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0 
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0 


1.0000 

0.1207 


1.02562 


0.000977 

0.1206 


3.023 

-4 

0.9993 

0.1207 

0.6204 

1.02483 

0.06727 

0.00293 

0.1205 


9.036 

-4 

0.9980 

0.1207 

0.6204 

1.02328 

0.06890 

0.01074 

0.1199 


3.265 

-3 

0.9926 

0.1208 

0.6205 

1.01742 

0.07563 

0.04199 

0.1177 


1.209 

-2 

0.9710 

0.1212 

0.6176 

0.99812 

0.10620 

0.1045 

0.1130 


2.736 

-2 

0.9282 

0.1217 

0.6038 

0.97058 

0.18106 

0.4795 

8.652 

-2 

7.814 

-2 

0.6972 

0.1241 

0.4834 

0.88725 

0.68001 

0.9795 

5.863 

-2 

9.047 

-2 

0.4660 

0.1258 

0.3361 

0.83598 

1.26512 

2.98 

1.084 

-2 

3.331 

-2 

8.412 -2 

0.1289 

6.445 -2 

0.76045 

3.19766 

8.98 

4.842 

-5 

2.433 

-4 

3.692 -4 

0.1311 

2.953 -4 

0.70911 

7.58946 

16.98 

2.886 

-8 

1.764 

-7 

2.187 -7 

0.1320 

1.775 -7 

0.69198 

11.8779 

24.98 

1.592 

-11 

1.061 

-10 

1.203 -10 

0.1323 

9.819 -11 

0.68540 

15.0246 

32.98 

8.487 • 

-15 

5.930 

-14 

6.406 -14 

0.1325 

5.243 -14 

0.68219 

17.3462 

40.98 

4.445 

-18 

3.194 

-17 

3.354 -17 

0.1325 

2.749 -17 

0.67994 

19.0513 

44.98 

0 


7.200 

-19 

7.492 -19 

0.0000 

5.003 -19 
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3.000977 
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1.0298 
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D. 00293 

1.772 -2 

1.506 -4 

0.9972 

1.777 -2 

0.9296 

1.0295 

0.07668 

3.01704 

1.760 -2 

5.428 -4 

0.9898 

1.778 -2 

0.9236 

1.0283 

0.08439 

3.04199 

1.709 -2 

1.989 -3 

0.9609 

1.778 -2 

0.8994 

1.0245 

0.1198 

3.1045 

1.611 -2 

4.423 -3 

0.9052 

1.780 -2 

0.8507 

1.0189 

0.2092 

3.4795 

1.126 -2 

1.175 -2 

0.6308 

1.785 -2 

0.5995 

1.0008 

0.8714 

3.9795 

6.948 -3 

1.266 -2 

0.3882 

1.790 -2 

0.3714 

0.9886 

1.750 

2.98 

9.810 -4 

3.785 -3 

5.460 -2 

1.797 -2 

5.281 -2 

0.9687 

5.158 

8.98 

2.593 -6 

1.797 -5 

1.439 -4 

1.802 -2 

1.403 -4 

0.9532 

15.069 

5.98 

9.117 -10 

8.289 -9 

5.005 -8 

1.804 -2 

4.941 -8 

0.9474 

28.054 

4.98 

3.156 -13 

3.314 -12 

1.749 -11 

1.804 -2 

1.712 -11 

0.9449 

40.929 

2.98 

1.084 -16 

1.253 -15 

6.008 -15 

1.804 -2 

5.884 -15 

0.9435 

53.740 

3.98 

3.709 -20 

4.596 -19 

2.055 -18 

1.805 -2 

2.013 -18 

0.9422 

66.505 

00 

as 

0 

8.724 -21 

3.794 -20 

0 
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0.99 
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0.99 

0.18816 
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0.86137 
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0.9 
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0.69386 
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0.5 
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0.5 
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0.5 

0.18816 
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32.98 
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8.98 
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0.87704 

463.12 

32.98 
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FIGURE CAPTIONS 


Fig. 1. 


Fig. 2. 
Fig. 3. 
Fig. 4. 

Fig. 5. 


Downward normalized interior radiance as a function of the cosine (y) 

of the zenith angle. Curves are given for various values of the optical 

depth within a very thick homogenous layer scattering according to the 

Rayleigh phase function. The upper curves are for £0^=0 . 99 and the 

lower curves for 0)^=0. 9. In each case the radiance has been multiplied 

by the factor (it / diffuse flux at depth r) , so that the variation 

of the radiance with y at different depths can conveniently be compared. 

These curves are for y o =0. 85332 (31.42°) and the incident plane 

containing the direction of the incoming beam ( = 0° for left half 

and (J) =180° for right half of figure). The solar horizon is at the 

left of the figure, the zenith is at the center, and the antisolar 

horizon is at the right of the figure. The incoming flux is normalized 

to unity across a plane at right angles to the incoming beam. 

Downward normalized interior radiance for to =0.5 and y = 0.85332. 

o o 

See caption to Fig. 1. 

Downward normalized interior radiance for w =0.1 and y =0.85332. See 

o o 

caption to Fig. 1. 

The ratio of the downward interior radiance at y = 1 to the downward 

interior radiance at y = 0.03785 (0 = 87.83°) as a function of the 

optical depth within the medium. Curves are given for u> q = 0.99, 0.9, 

0.5, and 0.1 and for y = 1, y = 0.85332 (0 = 31.42°), y = 0.53786 

o o o 

(0 = 57.46°), and y = 0.18816 (0 = 79.15°). 

O O o 

Downward normalized interior radiance for to = 0.5, y = 0.85332 

o o 

(0 = 31.42°), and y = 0.71392 (44.44°) as a function of the azimuthal 

angle <J>. 



Fig. 6. Downward normalized interior radiance for m = 0.1, y = 0.85332 

o o 

(0 = 31.42°), and jj = 0.71392 (44.44^) as a function of the azimuthal 

angle <p . 

Fig. 7. Upward normalized interior radiance as a function of the cosine (y) 

of the nadir angle. These curves are for y = 0.85332 (31.42°) and 

o 

the incident plane - 0° and 180°). See caption to Fig. 1. 

Fig. 8. Upward normalized interior radiance for y = 0.85332 and oj =0.5. 

o o 

See caption to Fig. 1. 

Fig. 9. Upward normalized interior radiance for y = 0.85332 and w =0.1. 

o o 

See caption to Fig. 1. 

Fig. 10. Upward normalized interior radiance for = 0.5, y Q = 0.85332 

(0 = 31.42°), and y = 0.71392 (44.44°) as a function of the azimuthal 
angle <J>. 

Fig. 11. The interior downward radiance or flux as a function of the optical 

depth within a very thick homogenous layer for (U q = 0.1. The upper 

curve is the downward flux for y = 1 and the next curve is for 

o 

y Q = 0.18816 (79.15)°. The next curves in order give the downward 

radiance for y = 0.18816 (79.15°) and the following five values for 
o 

y : 1(0.00°); 0.92226 (22.74°); 0.81141 (35.77°); 0.47254 (61.80°); 

0.03785 (87.83°). The azimuthal angle <)> = 0°. The dashed lines 
indicate an exponential fitted at large optical depths. 
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